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Abstract 

The 130 GeV gamma-ray hne based on tentative analyses on the Fermi-LAT data is hard to be 
understood with dark matter annihilation in the conventional framework of the MSSM. We point 
out that it can be nicely explained with two body decay of a scalar dark matter ((/)dm 77) by 
the dimension 6 operator suppressed with the mass of the grand unification scale (~ 10^^ GeV), 
C D \4>DM\^FpLvF'^'^ I^Q^r^, in which the scalar dark matter c^dm develops a TeV scale vacuum 
expectation value. We propose a viable model explaining the 130 GeV gamma-ray line. 
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A thermally produced weakly interacting massive particle (WIMP) would be the most 
promising dark matter (DM) candidate explaining 23 percent of the mass-energy density in 
the present universe It might be closely associated with a new physics at the electroweak 
(EW) energy scale beyond the standard model in particle physics. In this sense, the on- 
going indirect DM searches, which cover the EW energy scales, are expected to provide a 
hint toward the fundamental theory in particle physics as well as in cosmology. 

Recent tentative analyses 2I-Q] based on the data from the Fermi Large Area Tele- 
scope (Fermi- LAT) j^, ^ exhibited a sharp peak around 130 GeV in the gamma-ray spec- 
trum coming from near the galactic center (GC). The authors pointed out the gamma- 
ray excess could be a result from DM annihilation to a photon pair.^ Since DM should 
carry no electromagnetic charge,^ the annihilation, XX 11 is possible only through 
radiative processes. If the sharp peak of the gamma-ray around 130 GeV really origi- 
nates from DM annihilation, the annihilation cross section and the mass of DM would 
be {av)^^^^^ = (1.27 ± 0.32+°;^^) x 10-^7 cmVs (2.27 ± 0.57+°;^^ ^ 10"27 cmVs) and 
'^DM = 129.8 ± 2.41^3 GeV, respectively, when the Einasto (NEW) DM profile employed 
It is almost one order of magnitude smaller than the total cross section for the thermal 
production of DM needed for explaining the present DM density (~ 10~^ GeV/cm^), which 
is about 3 x 10~^^ cm^/s [l|. Indeed, the cross section of order 10~^^ cm/s is much larger 
than the expected estimation of one-loop suppressed processes, assuming a thermal relic DM 

On the other hand, at the tree level, DM may annihilate into other final states, e.g. 
W~^W~ ^ ZZ,bb,T~^T~ , , etc, which can produce secondary continuous 7-ray spectrum 
through hadronizations or final state radiations. Thus, one can derive the constraints on 
the DM annihilation cross sections for those channels from the Eermi-LAT 7-ray observation 
data. Current limits on those annihilation modes are at the level of O{10~^^ — 10"^^) cm^/s 



for mDM ~ 130 GeV. Eor more details, see Refs. llHlGj. Moreover, produced W and Z 



bosons can a. 
PAMELA 



so lead to a sizable primary contribution to the antiproton flux measured by 



17|, which provides another constraint of (9(10 ^^) cm'^/s on the DM annihilation 



into W+W- and ZZ 



18 



19| . Consequently, any annihilating DM model to explain the 130 
GeV 7-ray signal should also satisfy such limits. 

Actually, the full one-loop calculations of the neutralino annihilation into two photons 



20N22I show that the annihilation cross section of order 10~^^ cm'^/s is impossible in the 
region of 20 GeV - 4 TeV DM mass in the minimal supersymmetric standard model (MSSM). 
In fact, the neutralino in the MSSM can be annihilated also into one photon plus one Z 
boson through one-loop induced processes, and this gamma-ray can cause the excess of the 
observed flux. Unlike the case of xX ~^ II1 the emitted photon energy is estimated as 



^ In Rcfs. 0, 0|, it is argued that the gamma-ray line can be stiU explained with an astrophysical origin, 

associated with hard photons in the "Fermi bubble" regions. 

^ The possibility that DM is a milli-charged particle has been studied in Refs. |9l. llOf). 
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= mDM(l — ''^|/4?tidm) XX ~^ 1^- From the 130 GeV photon energy, hence, the DM 
mass of 144 GeV is predicted. [If the Z boson is replaced by an unknown heavier gauge field 
X, the DM mass can be raised more.] Even in such a case, the cross section of xx ~^ 1^ is 
just about 10-28 cmVs in the MSSM 



22 



23| . which is still smaller than 10 cm'^/s. 



Renouncing the possibility of thermal production of the neutralino DM required to explain 
the observed DM relic density, the annihilation cross sections {av)^-^^^^ and {av)y^y^^^z can 
be of order 10"^^ cm^/ s or even larger. In this case, however, the mass difference between the 



chargino and the neutralino should be of order 10 GeV or less |18l. |24| . Moreover, such large 
neutralino annihilation cross sections for the one-loop suppressed processes are necessarily 
accompanie d by much larger annihilation cross section into W^W~ of order 10"^^ cm'^/s or 



even larger 18|. As discussed above, the large annihilation cross section, (av) 



is 



12|-ll6|. 



constrained by the current Fermi-LAT limits on continuum photon spectrum 

If the neutralino is wino- or higgsino-like, one might think that the cross section of 
XX 11 could be enhanced by nonperturbative effects, called "Sommerfeld effect." It 
turns out, however, that the cross section {av)^^^^^ cannot reach 10"^'' cm^/s, unless the 
mass of the neutralino is of TeV or hundreds of GeV scales 25 



27l |. which is exceedingly 
heavier than 130 GeV. In addition, the nonperturbative effects on heavy wino- or higgsino- 
ike DM also enhance the annihilation cross sections into W^W~ and ZZ, {o'v)^^^w+w~ /zz 
ch are inevitably constrained by the Fermi-LAT continuum photon limits 



7 ^^vJ I 



19|. 



- |l6| and the PAMELA antiproton flux limits 
Thus, the 130 GeV gamma-ray is quite hard to explain with DM annihilation, if the 
framework is restricted within the MSSM: we need to consider a possibility of the presence 
of a new DM sector, introducing a new DM and its interactions with ordinary charged 
particles. The basic reason for the difficulty is the charged superparticles' masses circulating 
in the loop cannot be light enough to enhance the cross section, because they should be 
heavier than the neutralino DM. Hence, if a new interaction coupling between a new DM 
and charged particles is introduced, which is larger enough than the weak coupling, we may 
obtain the required cross section (av)^^^^^ with relatively heavy (130 -|- a fewxlO GeV) 
charged particles in the loop. Of course, the out-going interaction of the photons should be 
still given by the electromagnetic interaction. In order to reconcile the difference between the 
demanded cross sections for 130 GeV gamma-ray by DM annihilation and for the thermal 
relic DM, one may introduce two quite different interactions such that a photon annihilation 



interaction with (av) 



2 X 10 ^^cm'^/s is separated from the interaction explaining 



30|. However, we will not pursue such 



the thermal relic with ^{crv) ~ 3 x 10~2^cm^/s 
an ambitious job in this Letter. 

Instead, in this Letter we will discuss the possibility that the gamma-ray line at 130 GeV 
is explained by DM decay. By comparing the differential photon flux by DM decay ($dec) 
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with that by annihilation ($ann) 2m l22| 



dE^dVl 47r ' \2mDM/ Ji.o.s. ''^0 V 



V"^DM/ Ao.s. '^O V 



Phalo(^)' 






/ dE^ ' 






Po 


/ dE^ 



(1) 



one can estimate the decay rate F, needed for explaining the gamma-ray excess, where 
dN^eciann) / dE^ is the differential photon energy spectrum, pha\o{r) is the DM halo density 
profile, Pq ~ 0.4 GeV cm~^ is the local DM halo density, Tq ^ 8.5 kpc is the distance from the 
GC to the Sun and j^^^ ds is the integral along the line of sight (l.o.s.). The morphology of 
the signal from DM decay is linearly proportional to the DM density profile, while that from 
DM annihilation has the density square dependence. Consequently, the decay case tends 
to show a less steep increase of the signal towards the GC compared with the annihilation 
case, although the morphology of the signal still has uncertainty by the DM halo density 
profile itself. In addition, for decaying DM more 7-ray fiux is generically expected from the 
galactic halo compared to annihilating DM. However, the Fermi collaboration has observed 
no 7-ray excess from the galactic halo: it just reported the lower limits on the partial DM 
lifetime r^y jsl. Although the best-fit values for the lifetime are in tension with the limit 
from the Fermi collaboration, however, the required lifetime to explain the 130 GeV 7-ray 



signal marginally satisfies the experimental limit allowing 2a level error bars [33|. Moreover, 
there still exist the large uncertainty of the DM distribution around the GC, and also large 
statistical and systematic uncertainties at the moment. To confirm which scenario explains 
the 130 GeV 7-ray line, more improvement in observation is therefore essential in the near 
future. In Ref. [3^], it was shown that both of decaying and annihilating DM explanations 
similarly give good x^-^is for DM halo profiles more enhanced around the GC (with a > 1), 
compared to the original form of NFW profile (with a = 1). 

In Eq. ([T]), we set the DM mass in the decay case as two times heavier than the annihilation 
case to obtain the same resulting gamma-ray fiuxes around E^ = 130 GeV. Thus, the 
decay rate leading to the same gamma-ray fiux by the annihilation with (av)^^^^^ ~ 2 x 
10~^^ cm^/s is estimated as js], ^ 

F^^,, ~ 10-2^ s-\ (2) 

by which the life time of DM becomes sufficiently longer than the age of the universe (~ 10^^ 
s). Note that the required annihilation cross section for xx ~^ ^1 is approximately twice of 
that for XX ~^ 77! the required decay rate for x ~^ ^1 is also approximately twice of that 
for X 77- 

If the gamma-ray excess should be explained by DM decay, the sharp peak of the gamma- 
ray would imply two body decay of the DM, since the three body decay would make the 
spectrum broad and the intensity much weaker. In the case of x ~^ 7^, the emitted photon 
energy is estimated as E^ = {1 — m\/m^yi)m^yi/2. For E^ ^ 130 GeV, thus, the required 
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DM mass is around ttidm ~ 288 (1138) GeV for X = Z (1000 GeV), which is heavier than 
that in the annihilation case. Thus, in the decaying DM case the DM mass can be much 
heavier than 260 GeV, say upto C(TeV). 



In Ref. 35|, radiative DM decays to gamma-ray were extensively studied. For fermionic 
DM decay, the author considered the following renormalizable interactions: 

+ V^DM + V^Pi^ l^ + N [v^Pl + v^Pr] /S + h.c, 

where "(^''s and "?/"s denote the coupling constants, and Pl,r the projection operators. 
and S are superheavy vector and scalar fields with the masses my and m^, respectively, 
which radiatively mediate DM decay. and I indicate neutral and charged fermions, re- 
spectively. We suppose ttidm < 2m/ to disallow the three body decays of DM kinematically 
at tree level. Note that for producing the photons radiatively, the vector field (and also 
S) should carry an electromagnetic charge like the "X" or "y" gauge bosons in the SU(5) 



grand unified theory (GUT' 
for the fermionic DM ^dm 



. The interactions of Eq. ([3]) yield the estimation of the life time 
33: 



27 ^ /260 GeVy /^rtv_y ... 

" [E^{V9k4 - 9HW [ m..M J llO-GevJ ' 

provided my <^ m^. Thus, the life time can be of order 10^^ sec. e.g. for a charged gauge 
boson slightly heavier than 10^^ GeV. In many GUT models, however, the mass of the heavy 
gauge bosons carrying electromagnetic charges should be well above a few times 10^^ GeV 



for longevity of the proton [36[. Note that in the decaying DM CcLSG, clS discussed above, the 
DM mass can be much heavier than 260 GeV. In Eq. (jl]), thus, my could be raised upto 
10^^ GeV scale for fermionic DM of TeV scale mass. However, it is not yet enough to reach 
the conventional SUSY GUT scale (^ 2 x 10^^ GeV). On the other hand, if my S> m^, the 
life time of the fermionic DM becomes jssl 

in28 /260 GeVy / mg \4 

r,„,^,, " " ' lE^iy^y', - WH)]' V^^J I lO^^J • 

Thus, S with 10^"^ GeV mass still affects the gauge coupling unification, unless S composes 
an SU(5) multiplet with other fields. 

One could also explore the possibility of a scalar DM decaying to two gammas. In that 
case, however, the resulting effective dimension five operator {(p-QyiF^yF^" / M^) should be 
extremely suppressed for its life time of 10^^ s {ssj: the suppression factor (1/M^,) should be 
much smaller than 1/Mp, where Mp denotes the reduced Planck mass (~ 2.4 x 10^® GeV). 
Thus, let us consider the case that a scalar DM decays to two photons via the dimension 
six operator: 

-C., = c.,^^^F,^F^^, (6) 
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where the scalar DM, 0dm is assumed to develop a vacuum expectation value (VEV). With 
Eq. the decay rate of 0dm is estimated as 



0.85 X 10" 




(7) 

J V260 GeVJ 



Now we attempt to achieve this DM decay rate by constructing a simple supersymmetric 
(SUSY) model. The scalar DM, 0dm can be regarded as the scalar component of a superfield 
$. In this Letter, we simply assume that 0dm and also its fermionic superpartner $dm are 
non-thermally produced. Let us consider the following superpotential: 

W D /^iV$^ (8) 

where and $ are MSSM singlets, and k is a dimensionless coupling constant. The R and 
hyper charges of the relevant superfields are presented in Table [H We suppose that waterfall 
fields, which are not explicitly specified here, decay eventually to the scalar and fermionic 
components of $ through A^: A^"^ -> $dm0dm, N* ^dm'^'dm, and A^* 0dm0dm by 
the above k coupling and the corresponding "A-term." By the soft SUSY breaking yl-term 
corresponding to the k term in Eq. ([8]) and the soft mass terms in the scalar potential, 0dm 
and A^ can develop VEVs at the minimum: 

(0dm) ~ (N) ~ ^ ~ C(l) TeV. (9) 

K 

We assume that the mass of 0dm determined by the soft terms is about 260 GeV. The 
y4-term and the VEVs of 0dm and A^ (and also instanton effects) break U{1)r completely. 
However, the above singlet fields are still hard to be coupled to the ordinary MSSM fields 
carrying non-negative integer R charges: since A^ and $ carry positive fractional R charges, 
the ordinary R parity violating terms (and also the terms leading to the dimension 5 proton 
decay), which require the presence of a spurion field carrying the R charge of —1 (—2), 
should be extremely suppressed in this framework. As a result, such U{1)r breaking leaves 
intact e.g. the proton stability. 



Superfields 


N 




X 




Y 




U(l)ii 


4/3 


1/3 


1 


1 


2/3 


4/3 


U(l)y [=U(l)em] 








q 


-q 


-q 


q 



TABLE I: R and hyper charges of the superfields. The ordinary MSSM matter (Higgs) superfields 
including the Majorana neutrinos carry the unit (zero) R charges. 
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FIG. 1: Scalar dark matter decaying to two photons. By radiative mediations of the superheavy 
particles X and Y, the scalar dark matter component 0dm can decay into two photons with the 
decay rate of 0(10"^^) s~^. Similar diagrams contributed by the virtual superheavy scalar partners 
oi X, Y are also possible. 

0DM cannot be absolutely stable in this scenario, because of the following additional terms 
in the superpotential: 

W^^ = X^^^XY + MxXX' + MyYY", (10) 

where Mx and My are the mass parameters of the GUT scale (~ Mqut ~ 10^® GeV). 
The 1]{1)r and U(1)y charges of the superfields X^^^ and Y^'^^ are presented in Table [B If 
the model is embedded in an ordinary GUT such as SU(5) and SO(IO) GUTs, they should 
carry also SU(3)c and/or SU(2)i quantum numbers and be accompanied with other (colored) 
particles to compose proper irreducible representations and their conjugations of a GUT. On 
the other hand, if the gauge group is given by "flipped SU(5)" or just that of the standard 
model at the GUT scale, X^^^ and Y^^'' can still remain as singlet fields with q = ±1 without 



any other supplementary particles 37 1. 



The DM component, 0dm eventually decays into two photons with the desired decay rate. 
See Figure [H The loop-suppression factor, which gives a small Ces in Eq. (JTj), could be easily 
compensated by a relatively large VEV of 0dm (~ C'(l) TeV), the coupling A^^ of order 
unity, and a proper choice of charge q such that Eq. ([7]) can be fulfilled without any serious 
fine-tuning. In a similar way, 0dm can decay into 7Z: one photon is replaced with Z in the 
diagram of Figure [TJ For fn^^^^ = 260 GeV, the energy of this single photon is given by 

^ -!!M2M I ) ~114GeV. (11) 

2 \ mi 

\ <PDM / 

Considering the coupling difference and the reduced phase space factor, one can easily esti- 
mate relative decay rate: 

^4>uM^fZ ^Q26. (12) 
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Since only one photon is produced by the decay mode (pBU iZ, the final photon fiux 
is just 0.13 times of the 130 GeV gamma-ray fiux. Thus, we predict another peak around 

= 114 GeV although the 114 GeV second peak is much less significant than the primary 
peak around 130 GeV. 

On the other hand, the fermionic component of $, $dm cannot decay via a dimension 5 
operator: it is because $dm can decay to a photon plus a neutral fermion through a dimension 
5 operator only when a neutral fermion lighter than <I>dm exists as seen in Eqs. (jl]) and (|5]). 
However, we do not have such a neutral fermion. Note that in general a fermionic DM 
cannot decay to only two photons via a dimension 6 operator. 

Only with the interactions in Eq. ( ITOll and the above field contents, other terms yielding 
a dimension 6 operator, e.g. ^^LH^ cannot be generated from the renormalizable super- 
potential. By introducing more fields, however, ^^LH^ might be induced, since the given 
symmetries admit it: its suppression mass parameter, which is determined by the masses 
of such additional fields, would be model-dependent. Even in that case, the expected sig- 
nal — 7- l~H^ (or vH^) is kinematically constrained by the mass of the Higgsino: if 
the Higgsino is heavier than 0dm ? 0dm can decay only to three or more particles, which is 
much suppressed if the suppression factor is of order Mq^t oi heavier. The other possible 
decay channel 0ojy[ — )■ z/$dm might be detectable only in neutrino observatories such as 
Super-Kamiokande and IceCube. The current limits from the Super-Kamiokande data are 
tdm ~ 10^4-25 gg(, 3g|^ which is easily evaded if this decay mode originates from the dimen- 



sion 6 operator ^^LHu suppressed by at least MAut- In one year, the expected sensitivities 
of IceCube are around tdm ~ 10^^ sec. 38|, |39|. Consequently, it is difficult to observe 



the neutrino signal from (p^j^ — )■ i^$dm in near future. Of course, the bare superpotential 
can always contain ^^LHu/Mp. However, the Mp suppression results in too weak expected 
signals. 

If X and Y are accompanied with some other colored particles to be embedded in certain 
GUT multiplets, 0dm can decay into two gluons, gg, through the similar diagram to that 
of Figure [U two photons are replaced with two gluons. Then, these gluons can lead to a 
sizable primary contribution to the antiproton fiux. Thus, this decay mode 0dm — ^ gg is 
constrained by the PAMELA antiproton fiux data [l3]. Even though the direct constraint 
on gg channel has not been studied yet and the antiproton fiux constraints depend on the 
propagation models, we can estimate the limit as FV^ ~ ]^q26-27 g {^q^o. the limits on 
the W^W~, ZZ, hh, and qq channels [40 1. 

In conclusion, we have confirmed that the recently noticed 130 GeV gamma-ray line based 
on the Fermi-LAT data is hard to explain with DM annihilations in the conventional MSSM 
framework. We raised the possibility that it originates from two body decay of a scalar DM 
(0DM — ^ 77) by a dimension 6 operator suppressed with the mass of the grand unification 
scale {C D |0DMpi^/.^.i^^7M^uT)- The scalar DM needs to develop a VEV of TeV scale. We 
proposed a model realizing the possibility, in which superheavy particles with GUT scale 
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masses radiatively mediate the DM decay. 
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